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Abstract 

This article outlines the relative orbit control (guidance algorithm and its preliminary performances tests 
evaluation) that will be tested by the CNES Team on FFIORD (Formation Flying In Orbit Ranging 
Demonstration) onboard PRISMA mission. 

After a brief summary of the PRISMA mission context, the paper provides a full description of the rendezvous 
function involved in the approaching guidance experiment. This FFIORD onboard function is detailed in terms 
of on-board algorithmic method (basic algorithm and enhanced alternative), sensibility analysis used to 
construct manoeuvre plans and preliminary tests results. 


1 PRISMA MISSION AND FFIORD EXPERIMENT 

PRISMA is a “technology in-orbit testbed mission” [Ref. 1] for demonstrating formation flying and rendezvous technologies, 
as well as flight testing of sensor and actuator technologies. The mission is funded by the Swedish National Space Board 
(SNSB), and is undertaken as a multilateral project with additional contributions from CNES(Centre National d’Etudes 
Spatiales), DLR (Deutsche Zentrum fur Luft- und Raumfahrt) and DTU(Danmarks Tekniske Universitet). The prime 
contractor is the Swedish Space Corporation (SSC), responsible for design, integration and operation of the space and ground 
segments, as well as implementation of in-orbit experiments involving Autonomous Formation Flying, Homing, Rendezvous, 
and Three Dimensional Proximity Operations including Final Approach and Recede. 

Launch is expected in the first half of 2009 as secondary payload into a sun-synchronous LEO orbit of about 700 km, and the 
expected mission lifetime will be approximately 8 months. All flight operations will be controlled by SSC via the Kiruna 
ground station. This mission will consist of two spacecraft. One Main spacecraft of 140 kg with full 3-axis reaction wheel 
based attitude control and 3-axis delta-V capability and one simplified Target of 40 kg with coarse 3-axis attitude control based 
on magnetometers, sun sensors and magnetic torques. 

The FFIORD experiment [Ref. 3] led by CNES includes the RF sensor-based relative navigation and guidance functions, 
including proximity operations (close station keeping and translations), stand-by on relative “football” orbit, rendezvous and 
collision avoidance (This last function is developed by EADS Astrium within a partnership with CNES - [Ref. 5]). The 
onboard software is implemented as a Matlab/Simulink® library delivered to SSC and then autocoded into C code with the 
Real Time Workshop Embedded coder®. 

One of the objectives of the FFIORD experiment is to validate in flight the guidance functions. These functions are used 
mainly for closed loop manoeuvres, where they provide: 

• the desired Delta-V and the time of its execution for sparse manoeuvres, 

• the calculated reference trajectory in the case of continuous thrusting. 

In particular this experiment will allow evaluating the in-flight performances and autonomy level of these algorithms in 
preparation of future formation flying operational applications such as SIMBOL-X (a French/Italian/German X-ray telescope 
mission). 


2 RENDEZVOUS EXPERIMENT 

The purpose of the approaching experiment is to demonstrate the capability to perform a safe rendezvous of the MAIN 
spacecraft from a distance of about 10 km from the TARGET spacecraft relying on the RF navigation system and exercising a 
large degree of autonomy. The rendezvous will be achieved several times with different initial states and scenarios of 
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manoeuvres. For each test the performance will be evaluated based on 2 criteria: offset between the desired goal and the 
position reached, and deviation of the delta-V cost with respect to the expected value. 



The rendezvous must span over several orbits in order to minimize the delta-v cost and thus to allow a few rehearsals 
throughout FFIORD experiment. The duration of a rendezvous will be in the range of 12 hours for an average cost of 20 cm/s. 

In this experiment we do not consider essential to implement a folly autonomous system that does not need any planning from 
the ground since future formation flying like SIMBOL-X will not use such a feature. As soon as the rendezvous is initiated, our 
objective is to demonstrate the capability to run the whole rendezvous without ground intervention and adjust the satellite 
trajectory to circumstances even though the on-board computing of manoeuvres follows some uploaded plan. 


3 RENDEZVOUS BASIC ALGORITHM 

Consequently, the rendezvous algorithm has been designed to work with the onboard estimated relative state and a pre-defined 
manoeuvre plan { t k , AV k } (with k <10). The general principle is the following: an initial manoeuvres plan is uploaded by TC 

and the on-board algorithm update autonomously this plan with respect to the current estimated relative state X (t ) with the 

objective to get to the target X(t) on time t and without changing the dates of manoeuvres. The model of the relative 
motion used by the rendezvous function relies on Clohessy-Wiltshire equations (also called Hill equations) : 


with : 


X(t) = A{t - t c ) X{t c ) + B(t - t c ) r+ £ ■ A(t - t k ). AV k 

k=l 


(i) 


• A and B : Clohessy-Wiltshire matrices in the LVLH (Local Vertical Local Horizontal) frame of the target 
spacecraft. B represents the drag effect, it is optionnal. 

• y: differential drag acceleration (regarded as constant during all the rendezvous phase), 

• AV t ='(0,0,0,AV tx ,AV t j,AV t J: manoeuvre n°k in the LVLH frame of the target spacecraft. 


• t c : current date 

• t k : date of each manoeuvres of the plan ( n < 10) 

• t : date associated to rendezvous point 


If we split the matrix A(t ) in block matrices like: A(t) = 
Ar t ='(AK„AK„AV t ) , equation (1) can be rewritten more simply as: 


4(0 4(0 

K 4(0 A 4 (t)y 


and we rename AV k in the following way 


= c ('.o 

k=l 


( 2 ) 


with: C(t,t c ) = X(t)-A(t — t c ).X(t c )-B(t — t c ).y and M k 


'4(t~t k ? 

,40-t k )/ 
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Naturally, the role of the rendezvous function is to compute at each time step the serial of AV k that enable to get X{t) at t 
with the minimum consumption of propellant. Thanks to the expression (2) of the Clohessy-Wiltshire model, and a 
quadratic criteria for the minimization of the consumption, the problem solved on-board by the rendezvous algorithm can be 
reduced into a standard least square problem as follow: 


Ys M k- AV k - C(t,t c ) ( M DV = C(t,f c ) 

min(£||AK,f) “1 >™n(|nrf) 

k = 1 


with 


M = {M l ,M 2 ,...,M n _ l ,M n )e M 63n (9t) and DV= t (AV 1 ,AV 2 ,...,AV n _ l ,AV n )e Af 3nl (9t). 


The solution of this standard least square problem is linear as following : 


AF = Q n .C(t,t c ) with Q„ 


< M " 1 

M + 


if n > 2 
if n = 2 
if n = 1 


(3) 


( Note that another standard criterion like min IzLII AV k \\J instead of the quadratic criterion would have lead to a standard 

linear programming problem whose simplex solution is linear as well. The simplex alternative has the advantage to compute 
the manoeuvre dates which make it more autonomous but it is more difficult to implement on-board because of its sensibility 
to the manoeuvres dates. It is anyway another option to be studied). 

Finally, equation (3) sum up the on-board rendezvous algorithm mechanism: at each time step the algorithm compute C{t, t c ) 
and multiply it with Q n to deduce the serial of oncoming AV k . The Q n can be pre-calculated on ground at the beginning of 

rendezvous phase (the problem of its invertibility for the computation is considered on §5 ). Q n is different for each step of 
the plan depending on the number of remaining manoeuvres. The diagram below sum up the different phases of the rendezvous 
function involving the different Q during the experiment: 
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This onboard rendezvous guidance function needs an initial ground manoeuvre plan (mainly, manoeuvres dates t k and 
associated pseudo inverse matrices Q n , , Q n _i , . . . , Q 2 , Qi . the ground computed AV k are only used to monitor the onboard 
computed AV k ). Today, tool to be used to build this ground manoeuvre plan and its associated matrices will be certainly tool , 
so called DRAGON (Determination des Rendez-vous par Algorithmes Generaux d’ Optimisation Numerique), used at CNES 
to optimize rendezvous between near-circular and near-coplanar orbits. This tool, resulting from a tight collaboration with 
KIAM (Keldysh Institute for Applied Mathematics), has been involved for several years in phasing manoeuvres calculation in 
the frame of ATV mission (Automated Transfer Vehicle) and ATV - ISS rendezvous. 
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4 ENHANCED RENDEZVOUS ALGORITHM 


As indicated above, in FFIORD experiment we use a linear model, so called Clohessy- Wiltshire model, as a closed-form 
solution to the differential solution of relative motion on near-circular orbit. To obtain an enhanced rendezvous algorithm 
while keeping a robust and simple mechanism such as a linear mechanism, we introduce and justify in this paragraph a 
polynomial change of variable on the Clohessy- Wiltshire model that improve the accuracy of the model while keeping its 
linear structure. 

4.1 UNSTABLE VARIABLE AND LIMITS OF THE LINEAR MODEL 

To do so, let’s come back to the origin of the Clohessy-Wiltshire model. The relative dynamics of the main spacecraft is 
studied in the LVLFI frame of target spacecraft orbiting the Earth on a circular orbit (pulsation <n, radius R) without any other 
perturbations. Assuming these hypotheses, we represent the relative motion by the following dynamical system: 

X = F{X) (4) 

In order to study the motion of the main spacecraft around the target we develop the dynamic model (4) in the vicinity of the 
origin. Applying the Taylor development of F around 0 at the order 2 (lawful since F is C' in 0 thanks to the hypotheses ) it 
gives : 


(4 X = AX + F 2 (V) + o(x| 3 ) 


(5) 


where A 
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Note that in this development, the second term F 2 (X) is made of homogeneous polynomials of degree 2 in X coordinates. In 
particular F 2 (X ) = o|x| ). 

To obtain a clearer representation of the dynamics we write (5) it the basis that give A in its Jordan canonical form. With the 
linear change of variable: Y = M 1 X we obtain equation (5) in the Jordan basis of A : 

(5)^7 = 57 + F 2 (7) + o(7| 3 ) (6) 
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and F 2 (Y) = M 1 F 2 ( MY) .Note that F 2 (Y) is also made of homogeneous polynomials of degree 2 in Y coordinates (direct 
consequence of F 2 definition) .In particular F 2 (Y) = °\ Y \ 2 )- 
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The on-board algorithm presented in §3 solves only the linear part of the dynamics Y = BY . The result is a linear model that 
represents the first order of the dynamics in the vicinity of the origin. In the Jordan basis it gives: 

1 t 0 0 0 0 

0 1 0 0 0 0 

0 0 e~ ita 0 0 0 

0 0 0 e‘“ 0 0 

0 0 0 0 e~ ia 0 

0 0 0 0 0 e‘“ 

As expected this solution is simple, robust and easily computed on-board the spacecraft. One can recognize here the well 

known Clohessy-Wiltshire model or Hill model (summed up by Y = BY ). It represents an approximation of the dynamical 
system (4) at the order 1. Indeed the difference R between the approximated system (7) and the original system (4) can be 

written at the first order as a vector of homogeneous polynomials of degree 2 which is an order 2 in |F| : 

R{Y)= R{y x ,y 2 ,y^y A ,y 5 ,y 6 ) = F 2 {Y)+ 0 ^)= o\Yf) 

Consequently the system (7) and its associated order 1 solution should be relevant as long as the approximation of order 1 is 

\y\ 

true. Simulations results show that it turns out to be relevant only when < 0.15%(i.e. in the case of FFIORD it is 

R t 

relevant in a ball of 10 km). This unexpected feature has its origin in the linear solution. Indeed the linear solution of (7) has an 
unstable variable y x (/) = V, 0 + V, 0 ■ t that increases rapidly with the time. The consequence is that the variable y l (t ) makes 

the rest R (y, , y 2 , y 3 ,y 4 ,y 5 , y (> ) increase drastically with time. So it makes the linear approximation be quickly 
unrepresentative of the real dynamics (4). 

4.2 NO LINEAR SOLUTION TO SOLVE THE DYNAMICS UNTIL THE SECOND ORDER 

In order to increase the validity domain of the on-board solution and meet PRISMA requirements the most simple idea is to 
look for onboard resolutions of system (4) that solve it until the second order, i.e. that solves the following system: 

Y = BY + F 2 (y) (8) 

The most simple (and classical) way to do so is to find a well suited change of variable that reduces system (8) to an integrable 
system. Since the second order F 2 (Y) is made of polynomials, we look for a polynomial change of variable. Since F 2 {Y) is 
homogeneous of degree 2 the following type of change of variable is sufficient: 

SR 6 

(9) where P(Y) is made ot homogeneous polynomials of degree 2. 

Y > Y + P(Y) 

Unfortunately, as we show here after, such a polynomial change of variable can not reduce system (8) to an integrable system. 
Indeed, let’s apply this change of variable to the system (6). We obtain: 

(4) <=> (/ + DP(Y ))F = B -(Y + P(Y )) + F 2 (Y + P(Y )) + o|f + P(Y I ) where DP(Y ) stand for the Jacobian matrix 
of P(Y). 

Besides, F 2 (F + P(y)) = F 2 (y)+ o |f| ) and o|f + P(f)| )= o|y| ) since F 2 (Y ) and P(Y) are homogeneous 
polynomials of degree 2. 


M~\x{t)=Y{t) = 
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So the system (6) becomes (/ + DP(y))Y = B ■ Y + B ■ P{y) + F 2 (f)+ o(| 7| )• By multiplying each side by 
C r-DP{Y )) we obtain : 

Y - DP 2 (7) ■ Y = B ■ Y + B ■ P{Y) + F 2 (y) - DP{y) ■ B ■ Y - DP{y) ■ B ■ P{y) - DP{y) ■ F 2 ( 7 )- + o(f| 3 ) 

Besides DP(y) = o(lb) since P(Y) is composed of homogeneous polynomial of degree 2. 

So, with the change of variable, system (4) take the final form: 

Y = B ■ Y + B ■ P(Y)~ DP(Y) ■ B ■ Y + F 2 (Y) + o\y\') 

This system is integrable (in the sense “linearly integrable”) until the second order if B ■ P(y)— DP(y)- B - Y + F 2 (Y)= () 
i.e. if we can find an homogeneous polynomial of degree 2 in H 2 (vectorial space generated by the homogeneous polynomial 
of degree 2) such that B ■ P{y)— DP(y)- B ■ Y = —F 2 (Y ) (10). 

Let’s introduce the operator L B such that L b (y) = B Y . L B induces an endomorphism ad L B on H 2 defined by 
ad L b (Q) = \Q,L b \ = DL b .Q — DQ.L b for any Q in IF .(In this definition one can recognize the Lie bracket in the 
notation [,]). 

With these definition condition (10) can be rewritten as follow: (l (l) <=> ad L b (p(y)) = —F 2 (y) . It means that condition 
(10) is verified as far as we can find an antecedent to — F 2 {y) in IF by the endomorphism ad L /; . One can show (It is too 
long to be done here and it is not the main point of this explanation) that ad L B is not inversible on IF_ (the reason comes 

from the eigenvalue 0 of B ). So no antecedent to — F 2 (y) can be found. The consequence is that it does not exist a change 
of variable of the type (9) to integrate the system (4) until the second order (Note that this is not incompatible with the 
possibility to obtain high-order analytical solutions in particular cases. As an example, high-order analytical solutions can be 
developed in the case of periodic relative trajectories [Ref. 8 and 9]). 

Consequently, to sum up this paragraph, the dynamic can not be solved until the second order with a polynomial change of 
variables. So, a simple polynomial resolution of the second order dynamics can not be implemented on-board to enlarge the 
domain of validity of the PRISMA control algorithm. 


4.3 ENHANCED LINEAR MODEL OF THE RELATIVE DYNAMICS 

Nevertheless the previous development enables to find ways to enhance the first order on board resolution of the dynamics. If 
condition (10) doesn’t have a solution, it exists however some homogeneous polynomial Q of degree 2 in FI 2 such that: 

ad L b (q{y))-F 2 (y)= R(y 2 ,y 3 ,y 4 ,y 5 ,y 6 ) (11) where R{y 2 ,y 3 ,y 4 ,y 5 , y 6 ) is independent with respect to the 

unstable variable yi. 

Indeed we show with Maple (the Maple sheets are available on request to the authors) that the homogeneous polynomial of 
degree 2: Q(y) = M 1 ^(M Y) verify condition (11) if T is chosen as the classical curvilinear change of variable (at the 
second order) defined as follow: 


T: 




f 2 2 A 

R ’ R ’ 2R ’ R ’ R ’ R 


Consequently by applying on-board the polynomial change of variable C\ defined as following: 


9T 




on 


7 > Y+Q{Y) 

the dynamics and then solving the linear system (7), the system that is solved is totally equivalent to the following system in Y: 
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(i + dq(y))y = by + bq(y) 


( 12 ) 


This system differs from system (4) only at the second order. Indeed by multiplying each side by (/ — DP(Y )) we show that 
system (12) is equivalent to the following system Y = B ■ Y + B • Q(y)— DQ(y) ■ B • Y + o|t| ) . Thanks to the 
endomorphism ad L B it can be summed up by the equation: Y = BY + ad L b {q{y))+ o|f| ). Since Q verifies 
condition (11), system (12) is finally rewritten: Y = B Y + F 2 (y) + R(y 2 ,y 3 , y 4 , y 5 , y 6 ) + o|y| ) (13). 

This time, the difference between the on-board computed system (system (13)) and system (4) is i?(z 9 ,z 3 ,z 4 ,z 5 ,z 6 ). It is 
still a rest of the second order since R e H but it is not dependent on the unstable variable y l (/) . So this time it ensures the 

\y\ 

expected domain of validity over the whole ball - — <1% that occurs usually when only stable variables are involved. 

R 

Consequently, by applying on-board the change of variable, the accuracy of the model has been enhanced without changing the 
rendezvous method. It improve the rendezvous function in accuracy and it enlarges its validity domain. 


5 SENSIBILITY OF RENDEZVOUS FUNCTION 

In this paragraph, we develop the analytical sensibility matrices of the rendezvous algorithm. As in previous paragraph, we 
develop here the FFIORD rendezvous algorithm in the Jordan base to have a clearer representation of the dynamics, to put on 
relief the impact of the unstable variable of the dynamics on the sensibility of the algorithm and to facilitate the calculation. 

5.1 ON-BOARD AV COMPUTATION IN THE JORDAN BASIS 

With the linear change of variable: Y = M 1 X we obtain indeed the Clohessy-Wiltshire model (1) in the Jordan basis as 
following: 


Y{t) = M~ l X{t) -M~ x .A{t - 1 0 ).M.M~ l X{t 0 )+Y J M~ l .A(t - t k AV k 

k = 1 

i-e- Y{t) = j{t-t 0 ).Y{t 0 )+f j d{t-t k ).M -'AV k (14) 

k = 1 


To simplify the calculation, we transform equation (14) in a symmetric equation with respect to t and t 0 , by multiplying (14) 
on the left by the coefficient j(t n — /) . where t n is the date of last manoeuvre. It gives the equivalent equation for the motion: 


J{t n -t).Y{t) = j{t n -040+ZA -O-M- 1 AV k 

k = 1 


(15) 


We follow the same process as in the third paragraph to deduce from model (15), the on-board computation of the Clohessy- 
Wiltshire model. We develop it in the Jordan basis: 

• First we isolate the variables (a V k . ) and we rewrite j(t n — t k ).M 1 in block matrices £ M 3 3 (c) : 


AV k = 


IK, 


and j(t n ~t k ).M 1 = 


^ J\ (^n h: ) ^2 (hi h )1 
(hi ~ hi ) J 4 (hi — hi ) 


Then we gather the factor of the variables (a V k ; ), 


in matrices N k : N luk = 


r JA-tk) 

J 4 (hi ~ hi )/ 
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Then equation (13) can be written as following: 


(15 -t).Y{t) = j{t„ -t 0 ).Y{t 0 )+ 

k = 1 

• And finally by gathering the A V k in a global variable D V k and the matrices N k in a global matrix M k we 
obtain the compact following expression for the dynamics: 

J{t n -l).Y(l) = j(t n -t 0 )Y{t 0 )+M ia m^\ (16) 

Where: DV kn = : € M 3{n _ k+l) ^) and M pJ = ( Np yN p Je M 63p {%). 

Wj 

Note here that in equation (16), in opposition to the variable Y which is written in the Jordan basis to separate clearly the 
unstable variable from the others, the variable DV is kept written in the LVLH frame to facilitate the interpretations on the 
sensibility matrices elements. 

Then, we follow again the same process as in paragraph 3 to deduce from model (16), the onboard computation of the AV in 
the Jordan basis at each time step of the algorithm. At each date t the rendezvous algorithm receive from the navigation the 
main spacecraft state, Y (/. ) ( t c is the date associated to current time). Considering this state as the initial state, the algorithm 

compute the serial of AV in the variable D V that ensure the rendezvous at t in Y(t ) ( / is the date associated to rendezvous 
target point). Naturally, as in paragraph 3, the solution is divided into 3 cases that depend on the number of remaining 
manoeuvers: for t n _, < / < t n 

P 1 c P 

if p<n-2 

= Q r l/(t,-t\r(t)-j(t„ -y).r(Q]| with q p = • m, if P = n - 1 < i n 

M n.p if P = n 

Note in equation (17) that at each time step t c such that t l < t ;: <t the rendezvous algorithm compute in the variable 
DV p n the whole serial of (A V k ) <k<n that ensure the rendezvous. Actually in this serial of AV only the first one AV p will be 
applied by the spacecraft propulsion system while the other (A V k ) p<k<H will be recomputed after t p before being applied. 
Consequently the equation of the AV applied is actually: 



We develop here the sensibility of the AV computation with respect to the navigation errors and the propulsion errors. 
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For practical purposes, the DV n is computed on-board every seconds between t l and t p . In order to lighten the notations, 
we will consider that t c = t p (considering t c = t p — St seconds would just add to the equation a negligible term like 

At c -t P )y 

Consequently equation (18) that gives the expression of the applied A V is changed in the following equation: 

AC, = s{Q p ){j(t , -<>(<„) J 

If we add in this equation the error of propulsion SAV p at t it gives: 

AF, = &r p +s(Q p llj(t , -(>{(,)] (19) 

In this equation Y{t ) is the main spacecraft state at t . According to equation (16) its expression is theoretically: 

k = 1 

If we add the navigation error SY (t k ) at t k , the propulsion error SAV k and we replace AV k by its expression (19) it gives: 

y(t p ) = H P )+j(t p -t 0 H, ) + X N p _ u , .S{Q k ).[j(t n - - J(t n - t k ){Y{t k ) + SY{t k )) + SA V k ] (20) 

k = 1 


Consequently the expression of AV p with respect to the errors of navigation and propulsion is obtained by inserting expression 
(20) in expression (19). It gives: 


AV p = S\V p +S\Q p M 
-s(Q p U(t„-t r 

k = 1 

[t n -tJSY(t p ) J 

-j{t n -t k ){Y{t k )+SY{t k ))+SAV k ]\ 



Note that equation (21) is consistent with the strategy: AV depends only on the first p — 1 manoeuvres. It is then sensible 
only on the first p — 1 errors of propulsion ( ( 3 A V k ) k<p ) and on the first p errors of navigation ( (SY(t k )) k<p )• 


The onboard navigation provide the rendezvous algorithm with the variables ( Y{t k ) + SY{t k )) k<p , which is the state vector of 
the main spacecraft in the local LVLH frame attached to the target, at the time of computation t k . The sensibility matrix of the 
DV j with respect to the navigation errors is thus given by the Jacobian matrix of the variable DV \ with respect to SY{t k ) . 
It gives the following block matrix of M 3;j 6n (9t) : 


dSY^j 


0 


5 


DV/SY 


dAV„ 

^dSYitJ 


0 




0 

dAV n 


where 


a^(Oj 


dAV p 

dSY(t k ) 


-S[Q\j[t n -t) if p = k 

< S\Q p ).J\t n - t p ).N P _ U .S(Q k )j{t n -t k ) if p<k 

0 if p>k 
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The propulsion is represented in the rendezvous algorithm by the variables (A V k + SA V k \ <p , which is the AV executed by 
the main spacecraft in the local framework attached to the target, at the time t k . The sensibility matrix of the D V ] n with 
respect to the propulsion errors is thus given by the Jacobian matrix of the variable DV \ with respect to SA V k . It gives the 
following block matrix of M 3n 3n (91 ) : 


' dAV l 
dSAV l 


0 


s 


DV/SAV 


dAK 

JSAV t 


0 


A 


0 

dAV„ 

dSAvJ 


where 


dAV p 

dSAV k 


1 if p = k 

s(Q p \A',-t p \N r - U -S(Q i ) if p<k 

0 if p > k 


5.3 SENSIBILITY OF THE TARGETED STATE 

We develop here the sensibility of the targeted state Y(t) with respect to the navigation errors and the propulsion errors. The 
purpose of this paragraph is to put on relief the impact of the different errors on the algorithm performances with respect to the 
targeted state. By developing equation (20) until F(?),t date associated to targeted rendezvous, we obtain the expression of 
Y(l) with respect to the errors of navigation and propulsion as follow: 


Y{t) = j{t-t 0 )Y{t a ) + j{t-t l )^N luk .S(Q k ){j{t n -t).Y{t)-j{t n -t k ).{Y{t k )+ SY(t k ))+ SAV k ] 


k=\ 


(22) 


By the same argument as the one used for the calculation of the AV sensibility matrix , the sensibility matrix of the targeted 
state Y(t) with respect to the navigation errors is also given by the Jacobian matrix of the variable Y(t ) with respect to 
8Y (t k ) . It gives the following block matrix of M 6n 6 (91 ) : 


) Y(t)!SY 



where 


iY(t) 

dSY{t, ) 






Once again the same argument as the one used for the calculation of the AV sensibility matrix, the sensibility matrix of the 
targeted state Y(t) with respect to the propulsion errors is also given by the Jacobian matrix of the variable Y (/) with respect 
to SAV k . It gives the following block matrix of M 6h 3 (91 ) : 


' 57(0 N 

dSAV l 


s 


Y(t)/SAV ~ 


af(d 


where 


dAYjt) 

dSAV k 




5.4 CONCLUSION ON THE SENSIBILITY ANALYSIS OF THE RENDEZVOUS FUNCTION 

The numerical computation of the sensibility matrices for each case of rendezvous will obviously give the exact sensibility 
factor of each A V k and of Y(t) with respect to each error of propulsion or navigation. It will give in each case the 
interpretation of the behaviour of the rendezvous scenario with respect to the errors. 
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From an analytical point of view, we can already highlight two conditions on the manoeuvres dates t k that ensure 

homogeneous terms in the different sensibility matrices. Consequently these conditions will be sufficient to ensure a robust 
behaviour of the rendezvous algorithm. 

Indeed, all the coefficients of the sensibility matrices are products of terms like S{Q k ), j(t t — t m ) and N . 


On one side the terms of type J (/. — t m ) and N are of the same type since N 


p,i 


jStv-ta) 


. They contain 


exclusively terms of type constant or of type oij, — t m ) or of type e ' . These terms are homogeneous as long as the 

only divergent type which is ttif, — t ) stay close to 1 . Basically, this means that the strategy must ensure that there is not too 

many time (more than several orbits of the target spacecraft) between two manoeuvres. This give us the first condition for a 
robust behaviour of the rendezvous algorithm: 


CONDITION 1: the dates of manoeuvres t , must not be separated by more than a few periods of the target spacecraft orbits . 

The physical interpretation of this condition is that when you let too many time between two manoeuvres the error due to the 
navigation or the propulsion is integrated in the dynamics which contains a divergent variable. The risk is that the errors 
diverge with this variable and get very high. As the sensibility matrices suggest it, by avoiding too much time between the 
manoeuvres you ensure a regular control of the errors and a small dependence of your algorithm with respect to divergent 
error. 


On the other side the terms of type S (Q k ) are of another type. If they derive as well from the type ,/(/ ,. — t ) , they contain 
above all an inverse procedure that need some restrictions to stay homogeneous. The restriction is that the inverted matrices 

involved in the calculation of S (Q k ) must be indeed invertible. Two cases appear: if k < n — 2 Q k —M n k {m n k ! M n k ) 
and the matrix that should be invertible is M n k M k , while if k — n — 2 Q k = M n k 1 and the matrix that should be 
invertible is M k . In both cases a sufficient condition for the invertibility of the matrices is that the columns of M n k form a 
free family of vectors. Since M n k — {N n k ••• N n n ) a sufficient condition is that the columns of the (/V n k ) form a free 


family of vectors. Considering the expression of N n k an obvious sufficient condition is that (t n — t m ) i 1 {t n — t ) 


(t -t) 

n 

X n q / 



for 


m, q €E [l ,n\ (it ensures that the exponentials of different vectors are free with respect to each other). This gives us the second 
sufficient condition of robustness for the rendezvous algorithm: 


CONDITION 2: the dates of manoeuvres t , must not be separated by an integer number of half periods of the target 
spacecraft orbits. 

The physical interpretation of this condition is that locating two of the manoeuvres at an integer number of half periods make 
your algorithm loose some controllability parameters on the dynamics while its degrees of freedom (number of AV k x yorz ) is 

kept equal. Between the two manoeuvres, the control that is done on some of the dynamics variables is redundant. The 
rendezvous algorithm loose efficiency (loss of controllability parameters) gets instable (the dynamical system is no more 
invertible and the sensibility matrix raise drastically) and sometime even loose the total control of the dynamics (in the case of 
the last two manoeuvres at an integer number of half periods the system has no solution and at least one of the dynamics 
variable is out of control). 

An important point is that the two conditions listed before are sufficient conditions for a robust behaviour of the rendezvous 
algorithm. But they are as well necessary conditions in the case of a two manoeuvre transfer. Consequently, for any transfer, 
the last two manoeuvres should absolutely verify these 2 conditions. If not, a variable in the spacecraft dynamics may remain 
out of control at the end of the transfer. 
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6 PRELIMINARY PERFORMANCES RESULTS 

Monte Carlo simulations are currently performed to assess accuracy and robustness of the onboard rendezvous function. In 
order to lead these tests, an accurate simulator was developed by CNES under Matlab/Simulink®. This testbed environment 
makes possible to analyse the impact of the following perturbations : 

• perturbations such as geopotential’s higher orders, atmospheric drag, luni-solar attraction and Sun radiation pressure, 

• eccentricity of the TARGET: an eccentricity up to 0.004 is expected instead of a null eccentricity, as assumed in the 
Clohessy-Wiltshire equations (In simulation tests, eccentricity interval variation is equal to [0, 0.01]) , 

• navigation : the FFIORD navigation function with a statistical model of navigation errors is used inside this testbed. 
The FFRF subsystem is specified to provide the following information every seconds : inter-satellite distance 
(accuracy = 1 cm) and azimuth and elevation of line-of-sight between two satellites (specified accuracy < 1°). 
Considering these inputs, the navigation function delivers an estimated distance (X LV lh) with a residual noise lower 
than 2 mm for short and long range and an error of position on Y| V i 1 1- Zlvlh components of about 5 cm for short 
range (0 to 1000 m) and 60 cm for long range (1000 m to 10 000 m) [Ref. 3] , 

• Manoeuvre realization errors: MIB (Minimum Impulse Bit), splitting of AV on pulses, and so on. 


In this paragraph, we present preliminary results obtained on some rendezvous. : 

• Four manoeuvres Rendezvous [-200 m, 0 m, 0 m, 0 m/s, 0 m/s, 0 m/s] LV LH — > [-300 m, 0 m, 0 m, 0 m/s, 0 m/s, 0 
m/s] lvlh : 

o Position and velocity accuracy wrt theoretical TARGET : 


Standart Deviation on X =0.0091315 n 
Standart Deviation on Z =0.01301 1 m 


Mean on X=-299.9944 m 
Mean on Z=-0.0 10055 m 


Target position point cloud in the XZ plane 
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mean 


"M: 


x iPSi ’ 


0.03 - 
0.02 - 
0.01 - 
0 - 

- 0.01 - 
- 0.02 - 
-0.03 - 
-0.04 - 
-0.05 - 


Standart Deviation on X =0.0091 315 m Mean on X=-299.9944 m 
Standart Deviation on Y =0.016047 m Mean on Y=-0.0194 m 


Target position point cloud in the XY plane 
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0.04 

- : : :■ : : h 
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Standart Deviation on Vx =0.0049814 m/s Mean on Vx=-0.0049024 m/s 
Standart Deviation on Vz =0.0076645 m/s Mean on Vz=0.019645 m/s 

Target velocity point cloud in the VxVz plane 
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Standart Deviation on Vx =0.0049814 m/s Mean on Vx=-0.0049024 m/s 
Standart Deviation on Vy =0.0027802 m/s Mean on Vy=0.0032888 m/s 

Target velocity point cloud in the VxVy plane 
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o Maneuvers updated by on-board function: 

Maneuvers dates: t 0 + 30 s, t 0 + 6676 s, t 0 + 7436 s, t 0 + 8176 s and date associated to TARGET point : t 0 + 8206 s (Orbital 
period of the TARGET : 5920 s) 


AV 

components 

(LVLH) 

AVi (m/s) 

AV 2 (m/s) 

AV 3 (m/s) 

AV 4 (m/s) 

Mean 

St. dev 

Initial 

ground 

value 

Mean 

St. dev 

Initial 

ground 

value 

Mean 

St. dev 

Initial 

ground 

value 

Mean 

St. dev 

Initial 

ground 

value 

X 

0.006 

0.0006 

0.0056 

-0.002 

0.0035 

-.0058 

0.001 

0.0006 

4.5 10‘ 5 

-0.01 

0.0077 

0.00016 

Y 

-5.0 10‘ 5 

2.6 10‘ 6 

0.0 

-8.1 10‘ 5 

3.9 10‘ 5 

0.0 

0.0002 

0.0003 

0.0 

0.003 

0.0029 

0.0 

Z 

0.003 

4.1 10' 5 

0.0029 

0.009 

0.006 

0.0041 

0.008 

0.003 

0.0018 

0.023 

0.013 

-0.0028 


• Six manoeuvres Rendezvous [5000 m, 0 m, 0 m, 0 m/s, 0 m/s, 0 m/s] LV LH — »■ [1000 m, 0 m, 0 m, 0 m/s, 0 m/s, 0 
m/s] lvlh : 

o Position and velocity accuracy wrt theoretical TARGET : 


Standart Deviation on X =0.10423 m Mean on X=1000.1321 m 
Standart Deviation on Z =0.076918 m Mean on Z=-0.1 1992 m 


Target position point cloud in the XZ plane 
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Standart Deviation on X =0.10423 m Mean on X=1 000. 1321 m 
Standart Deviation on Y =0.028983 m Mean on Y=0.0042214 m 


Target position point cloud in the XY plane 
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Standart Deviation on Vx =0.0097541 m/s Mean on Vx=-0.0036467 m/s 
Standart Deviation on Vz =0.023717 m/s Mean on Vz=0.018642 m/s 


Target velocity point cloud in the VxVz plane 
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Standart Deviation on Vx =0.0097541 m/s Mean on Vx=-0.0036467 m/s 
Standart Deviation on Vy =0.0068194 m/s Mean on Vy=0.0049664 m/s 

Target velocity point cloud in the VxVy plane 
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o Maneuvers updated by on-board function: 

Maneuvers dates: t 0 + 30 s, t 0 + 18795 s, t 0 + 26201 s, t 0 + 33118 s, t 0 + 34104 s, t 0 + 34597 s and date associated to TARGET 
point : t 0 + 34627 s (Orbital period of the TARGET : 5920 s). 


Maneuvers cumulated mean values of AV l 4 : (0.15 m/s, 0.008 m/s, 0.25 m/s). 
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• Nine manoeuvres Rendezvous [-10 000 m, 0 m, 0 m, 0 m/s, 0 m/s, 0 m/s] LV LH — ► [- 100 m, 0 m, 0 m, 0 m/s, 0 m/s, 
0 m/s] L vlh : 

o Position and velocity accuracy wrt theoretical TARGET : 


Standart Deviation on X =0.092525 m Mean on X=-1 00.1 1 8 m Standart Deviation on X =0.092525 m Mean on X=-1 00.1 1 8 m 

Standart Deviation on Z =0.096292 m Mean on Z=-0.1 2607 m Standart Deviation on Y =0.01 883 m Mean on Y=-0.01 2353 m 

Target position point cloud in the XZ plane Target position point cloud in the XY plane 
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Standart Deviation on Vx =0.0082544 m/s Mean on Vx=0.0026658 m/s 
Standart Deviation on Vz =0.0082207 m/s Mean on Vz=-0.017966 m/s 


Target velocity point cloud in the VxVz plane 


Standart Deviation on Vx =0.0082544 m/s Mean on Vx=0.0026658 m/s 
Standart Deviation on Vy =0.0046365 m/s Mean on Vy=-0.0042728 m/s 

Target velocity point cloud in the VxVy plane 
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o Maneuvers updated by on-board function: 

Maneuvers dates: t 0 + 30 s, t 0 + 18138 s, t 0 + 33941 s, t 0 + 46452 s, t 0 + 54354 s, t 0 + 61926 s, t 0 + 62915 s, t 0 + 63902 s, t 0 + 
64550 s and date associated to TARGET point : t 0 + 64590 s (Orbital period of the TARGET : 5920 s). 

Maneuvers cumulated mean values of AV l 9 : (0.26 m/s, 0.009 m/s, 0.25 m/s). 


7 CONCLUSION 

This article has presented the rendezvous function involved in the FFIORD guidance experiment. The initial performances and 
robustness tests show a correct behaviour of the algorithm. They tend to show that the resulting accuracy on the targeted state 
reached by the algorithm is sufficient to comply FFIORD requirements . 

Nevertheless, this approach seems limited to near-circular orbits. And it is as well dependent on the ground segment for the 
computation of the maneuvers dates which limits the function autonomy. So, to prepare mission such as SIMBOF-X flying 
formation (high elliptic orbit, e ~ 0.8) other methods are studied, especially methods based on Model Predictive Control [Ref. 
6 and 4] method and state transition matrices (Yamanaka-Ankersen closed-form solution - [Ref. 7] ), to cope with the 
eccentricity effects and to investigate the onboard optimization of maneuvers dates . 


14/15 


8 REFERENCES 


[1] S. Berge, B. Jakobsson, P. Bodin, A. Edfors, and S. Persson, “Rendezvous and Formation Flying Experiments within the 
PRISMA In-orbit Testbed“, 6 th International ESA Conference on GNC, 2005 

[2] W. Clohessy and R. Wiltshire, “Terminal Guidance System for Satellite Rendezvous", Journal of the Aerospace Sciences, 
vol. 27, No 9, 1960, pp 653-658, 674 

[3] M. Delpech, P-Y. Guidotti, C. Mehlen and N. Perriault “CNES formation flying experiment on PRISMA mission” - 30 th 
Annual AAS Guidance and Control Conference - February 3-7, 2007 Breckenridge, Colorado. 

[4] R. Larsson, S. Berge, P. Bodin and U. Jonsson “Fuel Efficient Relative Orbit Control Strategies for Formation Flying and 
Rendezvous within PRISMA” - 29 th Annual AAS Guidance and Control Conference - February 4-8, 2006 Breckenridge, 
Colorado. 

[5] Ph. Cayeux, F. Raballand, J. Borde, J.-C. Berges and B. Meyssignac “Anti-collision function design and performances of 
the CNES formation flying experiment on the PRISMA mission” - 10 th International Symposium on Space Flight Dynamics - 
September 24-28, 2007 Annapolis, Maryland. 

[6] Tillerson M., Inalhan G., How J. P. “Co-ordination and control of distributed spacecraft systems using convex 
optimization techniques” - Int. Journal of Robust and Nonlinear Control 2002;12:207-242. 

[7] Yamanaka K., Ankersen F. “New state transition matrix for relative motion on an arbitrary elliptical orbit” - Journal of 
Guidance, Control, and Dynamics - Vol. 25, N° 1, January-February 2002. 

[8] Richardson D.L. and Mitchell J.W. “A third-order analytical solution for relative motion with a circular reference orbit” - 
Journal of the Astronautical Sciences - pp. 1-12 - Vol. 51, N° 1, January-March 2003. 

[9] Gomez G. and Marcote M.. “High-order analytical solutions of Hill’s equations ” - Celestial Mechanics and Dynamical 
Astronomy (2006) 94:197-211 


15/15 



